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Modelling nuclear effects in neutrino interactions in 1 GeV region 
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We evaluate nuclear effects in neutrino reactions in a framework based on a model proposed by Marteau with 
quasi-elastic and A production processes treated together. Nuclear effects include RPA corrections and A width 
modification in nuclear matter. 
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1. INTRODUCTION 

Description of A excitation region in neutrino- 
nucleus interactions is the most problematic in- 
gredient in Monte Carlo codes 1,. Experiments 
provide cross sections with a precision ~ 25% 
|2I,'{| . MC implementations are based on a combi- 
nation of Rein-Sehgal pion production model 4 
with Fermi gas model with different level of so- 
phistication in kinematical assumptions. The ba- 
sic dynamical rule is a factorization of interaction 
in two steps: (i) neutrino-nucleon interaction and 
(ii) re-interactions of outgoing particles inside nu- 
cleus. 

In this contribution we wish to present computa- 
tions done in a framework of more ambitious the- 
oretical scheme. Our model is based on Marteau 
model but it includes several modifications and 
also simplifications made in order to be able to 
compare better its predictions with experimental 
data. Marteau model describes on equal foot- 
ing A excitation and quasi-elastic processes. It is 
based on the non-relativistic Fermi gas with RPA 
corrections due to contact interaction terms and 
exchange of pions and p mesons [H] ■ Elementary 
2p — 2h excitations are also included. A mod- 
ification of A width in a nuclear matter is done 
using Oset's results ,7. ■ The model provides inclu- 
sive cross section for quasi-elastic and A excita- 
tion reactions and also contributions from several 
exclusive channels. The original model is rather 
complicated as local density effects are taken into 



account from the very beginning [5]. 
Our modifications are: 

i) We avoid complications due to local density 
profile of nucleus keeping a constant value of 
Fermi momentum kp — 225 MeV. The model 
becomes easier to handle and local density effects 
can be included at the very end in the MC ap- 
proach. 

ii) We adopt relativistic nucleons kinematics (we 
use relativistic generalization of the Lindhard 
function) . 

iii) We adopt original Oset results for the A 
width in nuclear matter to specific kinematics of 
neutrino-nucleus reaction. 

iv) We do not include 2p — 2h part as it requires 
further study UJ. 

2. FERMI GAS 

The basic cross section formula is: 
d 2 a 



dqdv 
where 



g| COSflgg 

32ttE 2 ^ 



(1) 



= 8(k^k v + k^k' v - 9f.uk ■ k' + ie^ a pk la kP)(2) 

is the leptonic tensor, E is the neutrino energy, 
k, k! denote lepton initial and final four-momenta, 
qH _ fcv — k'v — [y^ q) i s four-momentum transfer. 
The hadronic tensor is of the form 
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We use the hadronic weak current 
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Figure 1. CC Quasi-elastic cross section on 
free nucleons. Experimental points are taken 
from 



Figure 2. CC A ++ excitation cross section 
on free nucleons. Experimental points are taken 
from |2I3| and refer to v^p — ► ym + p. 
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where F± , F% , G a and Gp are the standard form- 
factors. In the frame q = (0,0, q) we calculate 
(x,y€{N,A}M- 



„oo / \ M x +M + v / M y + M + v 



2M X 



2M„ 



x(a° 0x (u,q)a° 0y {u,q)R c xy (v,q) ( 5 ) 
+$ x { V ,q)$ y {v,q)R l xv ( V , q y 



a° Jv,q)=F 1 (Q 2 )-F 2 (Q 2 ) 



2M(M X + M + v) 



(G) 



/3°o>,9) = 
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(7) 



M x + M + v 2M M x + M + v J ' 
etc. The free Fermi gas is characterized by 
Rna — Ran = 0, (8) 



Vol 



lmll%_ h (i>, q), 



(9) 



■ f^NA^Vol 



J&aM = -(i^) 2 — XmU° A _ h (u,q), (10) 



' f-xNN 7T 
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2k p 
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9 (2nfJ -"(s-Mlr + M A Tl 
T A = PBL ■ T nN - 22Tm(S A ) 



(11) 
(12) 

(13) 

(14) 
(15) 



where PBL € [0, 1] is the Pauli blocking factor 
defined as follows |7]: 



PB = 



\PA\\q cm \ - y/sEp + E A E 

C 

2 |p*A 1 1 <fcm I 



(16) 
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CC v Quasi-Elastic Cross Section on 16 in Fermi Gas Model CC v A Excitation and ji Production Cross Section on ,6 in Fermi Gas Model 




Figure 3. CC Quasi-elastic cross section on 
16 O in Fermi gas model. 



Figure 4. CC n production cross section on 
16 O in Fermi gas model is smaller then A excita- 
tion cross section since A can decay without 7r's 
in final state 



PBL = 1 if PB > 1, PBL = if PB < 0, other- 
wise PBL — PB. (Ea,pa) is A 4— momentum, 
(E cm ,q cm ) is pion (from A decay) 4— momentum 
in the center of mass frame, 



(v^) 3 Ma 



(17) 



T = 115MeV, M A = U32MeV. Im(E A ) de- 
scribes nuclear effects in the form of extra contri- 
butions to the A width from channels A — > irN, 
NA -> NN, NNA -> NNN. 



Im(E A ) = Jm(E^) + Jm(S 



NN,NNN\ 



(18) 



Iiti(Ea) was calculated by Oset-Salcedo |Z] in 
two kinematical situations: pion-nucleon and 
photon- nucleon scatterings. Because the kine- 
matical region for neutrino induced reaction is 
different (y — g 2 < 0) wc adopt an approximation 



Im(YjA) v — Im(E A ) 7 + (Zt7j(I]a)- 



(19) 



with Tm(Y>A)-y, Ito(Sa)tt taken from Oset- 
Salcedo. Our prescription introduces an increase 
with respect to the value of Im(EA) 7 by amount 
of 5-10%. 



Normalization factors are checked by perform- 
ing the limit kp — > in which we recover quasi- 
elastic (Fig. 1) and A excitation (Fig. 2) cross 
sections on free nucleons. In the case of A ex- 
citation Marteau model provides a sum of cross 
sections over isospin states i.e. a sum of and 
A + productions. We recover A ++ production 
cross section assuming isospin branching ratio 
rule cr(A ++ ) = 3er(A + ). We perform a com- 
parison for this particular channel of pion pro- 
duction since it is known that in the reaction 
v^p — > /j,~ptt + resonance contribution is domi- 
nant. Experimental points are taken from papers 
|10I11I12] in the case of quasi-elastic reaction and 
|2I3| in the case of A ++ production. 

In calculating H^" we adopted an approxima- 
tion in which calculation of sums over hadronic 
spins is done in the limit in which target nucleon 
is at rest. In the case of quasi-elastic one can 
verify by explicit computations that this approx- 
imation is very good (on the level of 2-3%). In 
the case of A excitation computations of response 
functions done by Marteau |Sj show that in the 
relevant kinematical region the approximation is 
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also very good. 

Results in the Fermi gas model are presented 
in Fig. 3 (quasi-elastic reaction) and Fig. 4 (A 
excitation and tt production). In these and next 
computations we have assumed that nucleus in 
question is le O i.e. A = 16. As mentioned in 
the introduction we do not take into account 16 O 
density profile keeping a constant value of Fermi 
momentum. 

3. RPA CORRECTIONS 

RPA equations read (separately for c,l,t; we 
omit the arguments v and q of all the functions 
below) : 



RPA Effects in CC v Quasi-Elastic Differential Cross Section on 16 0; E =1 GeV 
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Figure 5. CC quasi-elastic cross section on 16 O 
with RPA corrections. 



itaa = n A _, l +n A _, l VAArnAr A +n A _ ft VAAnAA, 

IIjvA = n jv _ /l VjvjvIIjvA + IIjv.^VjvAnAA, 



iiaat = n A _ /l VAjviijvjv + n A _ ft VA A n A jv- 



(20) 



Figure 6. CC differential quasi-elastic cross 
section on 16 with RPA corrections; neutrino 
energy is lGeV. 



The solutions are found to be 

n NN - n° N _ h (i - VaauI^d- 1 
u AA = u A _ h (i-v NN u%_ h )D- 1 

IIjva = IIaat = VN A U A _ h U° N _ h D^ 1 
where 

D = (1 - v NN n%_ h )(i - V AA U%_ h ) 







V 2 17° TT 

After substitution 



R 



i,t ,UnAs 2 VoI 
a a = -(7 ) — ZmR 

JttNN tt 



i,t 

AAi 



D c,l,t Vol c i t 

r nn = XtoII„' 

TT 



N N ' 



T l,t 



nl,t _ r>l,t _ IttNA V ol ^ 
n NA ~ n AN — t „ -L'mLL NA 

JttNN TT 



(21) 



(22) 

(23) 
(24) 
(25) 



we obtain the final expression for inclusive the 
cross section. 
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RPA Effects in CC v Quasi-elastic Differential Cross Section on 16 0; E =1GeV 
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Figure 7. CC differential quasi-elastic cross 
section on 16 O with RPA corrections in decom- 
position to charge, longitudinal and transverse 
parts; neutrino energy is lGeV. 
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H V z — g z — TO^ 



— v 2 + q 2 
' A 2 -m 2 - 
vA 2 - i/ 2 + q 2 . 



(27) 



(28) 




G 



If^f = 2, /' = 0.6, g' = 0.7, g" = 0.5, </' 

J ttNN 

0.5, A T = lOOOMeV, A p = 1500MeV. 



4. EXCLUSIVE CHANNELS 



RPA Effects in CC v„ n Production Differential Cross Section on lb O 



Transverse Part with RPA Effects 

Transverse Part in FG model - 
Longitudinal Part with RPA Effects - 
Longitudinal Part in FG model 



In identifying contributions to exclusive chan- 
nels we distinguish quasi-elasic processes and re- 
actions with and without pions in the final state. 
For this purpose we split A width and response 
functions in two parts according to 



RPA Effects in CC v n Production Differential Cross Section on ,B 




Figure 9. CC differential 7r production cross 
section on le O with RPA corrections; neutrino 
energy is lGeV. 
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Figure 10. CC differential tt production cross 
section on 16 O with RPA corrections in decompo- 
sition to longitudinal and transverse parts; neu- 
trino energy is lGeV. 



culated by making substitutions 



R 



i,t 

NA 



\ n N&.)ir 



f-xNN n 



■(Jmn^ A ) w 



R 



i,t 

AA 



(-RjVJv)t 



Vol 



(32) 



•(^) 2 -(I^ A )„ 
JnNN tt 



Xm(U NN ) n = —2 , (33) 



\D\ 



-.00 



PBL ■ T^n - 2Im(Y, A ) 



ImU° A _ h =Im(n A _ h ) 7T +Im(U A _ h ) N 



lm(Il A _ h ) 7I (i>,q) = 



16 Ml 



{s-M A y + M A T\ 



d 3 P 



(29) 
(30) 

(31) 



-(Im(U NA ) n + lm{W AN )^) 



(34) 



(V NA (Re(U%_ h ) - ^in^J 2 )^!!^,^ 
\D\ 2 

2m(n A A)7r = rjjp (35) 



5. RESULTS AND DISCUSSION 



Contributions to pion production come from all 
four ImH X y sectors of the theory. They are cal- 



Our results concerning relevance of RPA are 
summarized in a series of plots. 
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In Fig. 5 total cross sections for CC quasi- 
elastic reaction are presented. We compare three 
situations: Fermi gas without RPA correlations, 
RPA correlations without taking into account 
A — h excitation and full RPA computations. It is 
seen that RPA does not introduce much change in 
the total cross section. Inclusion of A — h excita- 
tions increases cross section for neutrino energies 
above 0.8Gey and one arrives at results very close 
to those for Fermi gas. 

In Fig. 6 CC quasi-elastic differential cross sec- 
tion is shown for neutrino energy E — IGeV. 

A well know result is reproduced: the quasi- 
elastic peak is cut and some increase of differ- 
ential cross section is observed at bigger values of 
energy transfer 

In Fig. 7 transverse, longitudinal and charge 
contributions (according to spin operators 
present in the transition amplitude) are dis- 
tinguished in the CC quasi-elastic differential 
cross section and patterns of RPA modifications 
for them is seen. The transverse contribution 
dominates and the longitudinal one is negligible. 

In Fig. 8 total cross sections for CC tt produc- 
tion are shown. We compare two situations: with 
and without RPA correlations. Inclusion of RPA 
correction reduces the cross section significantly 
by about 25%. 

In Fig. 9 CC tt production differential cross sec- 
tion ^7 is shown for neutrino energy E — IGeV. 

In Fig. 10 in the CC tt production differential 
cross section longitudinal and transverse contri- 
butions (according to spin operators present in 
the transition amplitude) are distinguished and 
patterns of RPA modifications can be followed. 
There is a substantial reduction of the dominant- 
ing transverse part at bigger values of energy 
transfer. The longitudinal part is in general re- 
duced (order of 50 %) but some increase is also 
seen for smaller values of energy transfer. 

To summarize: we have presented a scheme to 
calculate nuclear effects in neutrino-nucleus in- 
teractions. As a candidate to implement in MC 
codes the model presented in this paper has to 
be supplemented with a non-resonant contribu- 
tion to tt production, perhaps after It is 



unclear how much of nuclear effects described by 
RPA and A width modification are covered in 
FSI (Final State Interactions) models of existing 
MC codes. This has to be understood in order to 
avoid double counting. Finally, as mentioned in 
the introduction, local density effects can be in- 
vestigated by repeating present computations for 
several values of Fermi momentum and by taking 
an appropriate average. 
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